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Abstract 

Using the CLEO III detector at CESR we study e + e~ collisions in the center-of-mass energy close 
to, or above, A^A^ production threshold. We search for evidence of A^A^ resonance production 
and set upper limits based on inclusive hadron production as a barometer of A°A° production. 

PACS numbers: 14.65.Fy, 13.66.Bc, 13.30.Eg 
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I. INTRODUCTION 



The A°, consisting of 6, u and d quarks, is the lowest-lying 6-flavored baryon, about 
which comparatively little is known. Recently the CDF collaboration reported an improved 
measurement of the A° mass |l( of 5620.4 ± 1.6 ± 1.2 MeV. The lifetime has long been 
measured to be somewhat lower than theoretical expectations There is, however, no 
measurement available on the direct production of exclusive A°A° in e + e~ annihilation. Such 
events would be very useful for establishing absolute branching ratios and other properties. 
CLEO has accumulated data using e + e~ collisions in the center-of-mass energy range from 
11.227 to 11.383 GeV, close to or just above the A°A° production threshold. It is possible to 
observe a resonant signal, similar to the T(4S) for B + and B° mesons, or just an increase in 
relative production above threshold. We report here limits on such resonant or non-resonant 
production. 

II. DATA AND MONTE CARLO SIMULATED SAMPLE 

The CLEO III detector is described in detail elsewhere 0|- The inner part of the 
detector is surrounded by a 1.5 T solenoidal magnetic field. From the region near the e + e~ 
interaction vertex radially outward it consists of a silicon strip based vertex detector and 
a drift chamber used to measure the momenta of charged tracks based on their curvature. 
Beyond the drift chamber is a Ring Imaging Cherenkov Detector, RICH, used to identify 
charged hadrons, followed by an Electromagnetic Calorimeter, EC, consisting of nearly 8000 
Csl crystals. Next to the EC there is the solenoidal coil followed by an iron return path with 
wire chambers interspersed in 3 layers to provide muon identification. 

This study is based on the total 710 pb -1 data sample that was acquired at 3 MeV 
intervals between center-of-mass energies, Eqm, of 11.227 GeV to 11.383 GeV, to be close to 
or above threshold for A°A° production. The luminosity in each of these scan points varies 
from 14 to 20 pb _1 . In addition, there are data points taken at a Ecu of 11.150 and 11.203 
GeV, respectively. The two data points with lowest and highest energies have integrated 
luminosities of 70 and 120 pb _1 , respectively. We also use data taken in the four-flavor 
continuum below the T(45) to measure the bb cross section above the T(AS). 

For the Monte Carlo, MC, study of the high energy data, we generated five times more 
hadronic qq events than at each beam energy contained in our data sample. Events were 
generated separately for "light" four-flavor continuum (c, s, u, d) and bb continuum events 
and then combined in the expected 10:1 ratio absent any resonance production. The decay 
channels and the branching fractions of the A& are less well known than the B° and B + 
mesons. We list the A& decay modes and branching fractions we used for the signal Monte 
Carlo in Table |1J For the A° — > h^l^v branching fraction we re-scaled the B° — > Xlv 
branching fraction by the ratio of lifetimes, r(Aft)/r(B°). The entries denoted by *qq* indi- 
cate that the processes are generated using a fragmentation process for the quark-antiquark 
pair. 
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TABLE I: A b decay modes and branching fractions used in the Monte Carlo simulation. 



Decay modes Branching fraction (%) 



A b 


-> A+e~v e 


8.4 


A b 


-> A+fi-u^ 


8.4 


A 6 




4.2 


A b 




1.0 


A b 


-> A+al- 


2.1 


A b 


- A+D7 


2.1 


A b 


- A+Z?r 


4.2 


A b 


-> A?? c 


0.1 


A b 


-> AJ/V> 


0.5 


A b 


— > A+7r~7r + 7r~ 


2.1 


A b 


— > Ai^ 7T _ 7r _ 7r + 7r _ 


2.1 


A b 


-> p + D°ir- 


2.1 


A b 


— > A+ * tfu* 


44.9 


A b 


— > S+ * cfu* 


8.4 


A b 


— * 0+ * cfe* 


7.3 


A b 


— ► p + * cftl* 


1.1 


A b 


— > * tfu* 


1.0 



III. EVENT SELECTION 

The major backgrounds to A b are non-66 type hadronic events, two-photon events (e + e~ — > 
e + e~X) and t + t~ pairs. To suppress these backgrounds we require the following hadronic 
event selection criteria: 

(1) At least five charged tracks; a track candidate is acceptable if it is a cosine with respect 
to the beam line of less than 0.9 and has at least half of the potential tracking chamber hits 
along its length. This requirement rejects 81% of the t + t~ pairs. 

(2) The total visible energy, E vis , is required to be greater than the beam energy, E beam . 
E vis receives contributions from both charged tracks and unmatched neutral energy clus- 
ters greater than 30 MeV. This requirement helps suppresses two-photon events. Fig. G^a) 
shows the E vis / E beam distributions for data, five flavor Monte-Carlo continuum and simu- 
lated two-photon events 0. Imposing the requirement E vis > E beam reduces the two-photon 
background by 75% with a small (3%) loss of hadronic events. 

(3) The ratio of the 2 nd and th Fox- Wolfram moments, R%, is less than 0.25 jjj. Fig.[TJb) 
shows MC simulated distributions of R2 for both 66 and non-66 continuum events. Both 
areas are normalized to unity. Requiring R2 < 0.25 selects the more spherically shaped 
events in momentum space and greatly enhances the 66 fraction, by rejecting 65% of four- 
flavor continuum events while losing only 8% of the 66 events. 

To subtract four-flavor continuum background we use data taken at a Eqm 30 MeV below 
the T(4S) mass. Since we make a specific cut on R 2 we need to take into account that the 
shape of the R2 distribution can change when the Ecu changes. The R2 distribution from 
below- Y^S 1 ) data is compared with the distribution using data taken in the A b scan region 
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FIG. 1: (a) E V i S / Eb eam above A& threshold data (triangle), five flavor continuum MC (solid) and 
simulated two-photon events (circles), (b) R 2 distribution for 66 (dashed) and non-66 (solid) type 
events. 

in Fig. E(a). The data are normalized by luminosity and 1/s, where s is the square of the 
center-of-mass energy. The distributions differ in two respects. The first is the obvious 
enhancement at small R 2 values in the A& scan region giving evidence for bb production. The 
second is the disagreement in shape at values of R 2 > 0.5, where bb production is absent. 

We confirm this change in shape with energy by comparing T(4S) "on- resonance" data 
and below-T(lS') resonance data {Ecu =9. 43 GeV) in Fig. 12(b). The subtracted spectra 
show an anomalous peak near R 2 = 0.5. The number of events in this peak can be as large 
as ~30% of the total number of bb events in higher Eqm region. Thus, it is important 
to transform correctly the below-T(4S') resonance data in order to correctly subtract the 
background when we apply a tight R2 requirement. Simple kinematic considerations suggest 
that R' 2 (E') I 'R 2 (E) ~ E'/E, where E' > E. The boundary considerations that at R 2 
values of both and 1 the initial and corrected distributions be equal, result in a simple 
parameterization of the corrected, or "boosted" R 2 distribution: 

R' 2 (E') = ^R 2 (E) + (l " f ) Rl(E) • (1) 

This expression describes the energy dependence of the R 2 shape excellently. In Fig. |3] we 
compare the boosted R 2 distribution for below- Y^S") data, normalized by luminosity and 
with the same distribution for the high energy data. The distributions match above R 2 
of 0.5, as required. 

We have several strategies for observing the production of A°A° events. One possibility 
is to look for enhancements in the (1) bb cross-section. Another is to look for an increase 
in (2) A or (3) anti-proton production. We don't use protons because there is a large 
background rate from hadron interactions in the beam pipe and from residual beam gas 
collisions. A's are promising because we expect that A£ — ■> A C X has a large branching 
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FIG. 2: The R2 distribution above A& threshold compared with below-T(4S) data (a) and T(45) 
on resonance data compared with below-Y(lS') data (b). Circles show the subtracted distributions. 



ratio, ~96% and A+ — > AX is approximately 50%. Detecting anti-protons is very promising 
because A° decays always produce either one proton or neutron. In the case of non-resonant 
A°A° production we can expect that the cross-section will increase from zero at threshold 
to some constant fraction of the total bb cross-section. In order to ascertain an optimal 
search strategy, we assume this fraction is 7.9%, as predicted by the JETSET 7.3 Monte 
Carlo model |7j. This is consistent with the PDG value for bb — > baryon of 10% [8(. Further 
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FIG. 3: R2 Distribution at one energy point above A& threshold compared with below-T(45') after 
the boost (data). 
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support for this value comes from the ratio of A C A C to cc rates. As input to this estimate 
we use a measured value of i3(A+ — > pK~n + ) x cr(A+) = (10.0 ± 1.5 ± 1.5) pb from 
our below-T(45') continuum data sample. We take the cc cross section as 4/10 of the total 
hadronic cross section, implying <r(cc) = 1.12±0.02 nb and we use the PDG mean value 
for £(A+ -> pK-n + ) = (5.0 ±_1.3)% 0, yielding the ratio or AX/cc = (8.9 ± 3.0)%. 

The relative size of the A°A° component for our different search strategies is shown in 
Fig. Ufa). Here we normalized the MC simulated five-flavor visible hadronic cross section to 
unity, defined here as "continuum" udsc and b, and then added the signal A°A° to the total 
udscb cross section(i.e., the A°A° enhancement here represents an additional 7.9% above 
expected inclusive bb hadronic cross-section, rather than simply presenting an additional 
channel available to bb hadronization) . A's have the highest relative yield closely followed by 
anti-protons. We optimize our search criteria by maximizing signal divided by square root 
of the background, S/ VB, for our different search methods. The results are summarized in 
Fig-EJb), where we show the statistical significance for signal we obtain for different analysis 
strategies for different A°A° cross-sections (statistical errors only). 




bb p A bb p A 



FIG. 4: (a) Relative yield of the udsc (lower), b (middle) and A& (upper) visible cross section for 
the inclusive selection of bb , p and A assuming a 7.9% increase of the total bb cross section above 
A^A^ threshold, (b) Signal/VBackground for different analysis strategies and cross-sections. 

Our studies indicate that baryon production (namely anti-protons and A's) is the most 
sensitive measure of A°A° . However, the systematic uncertainties in A& — > protons and 
Afe — > A diminish their sensitivity relative to inclusive bb production. We also considered 
identifying A's and protons with an additional lepton in the event but these methods offer 
less significance. The efficiencies for detecting hadronic events, and more importantly, for 
detecting events with one or more protons are listed in Table [H] their evaluation will be 
discussed in more detail in the next section. 

We use both charged particle ionization loss in the drift chamber (dE/dx) and RICH 
information to identify anti-protons. The RICH is used for momenta larger than 1 GeV. 
Information on the angle of detected Cherenkov photons is translated into a likelihood of a 
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given photon being due to a particular particle. Contributions from all photons associated 
with a particular track are then summed to form an overall likelihood denoted as L{ for 
each particle hypothesis. To differentiate between kaon and proton candidates, we use the 
difference: — log(£^) + log(£ pr . oton ). This cut is set at -4. To utilize the dE/dx information 
we calculate ok as the difference between the expected ionization loss for a kaon and the 
measured loss divided by the measurement error. Similarly, a pro ton is defined in the same 
manner using the expected ionization for a proton. 

We use both the RICH and dE/dx to select anti-proton candidates in the following 
manner: (a) If neither the RICH nor dE/dx information is available, then the track is 
rejected, (b) If dE/dx is available and RICH is not then we insist that proton candidates 
have PIDdE = &\ — & proton < ( c ) ^ RICH information is available and dE/dx is not avail- 
able, then we require that PIDrjch = — log (£/<:) + \og(£ prot on) < —4. (d) If both dE/dx 
and RICH information are available, we require that (PIDdE + PIDrich) < —4. 

A candidates are formed from a pair of oppositely charged tracks one of which is consistent 
with a proton or anti-proton hypothesis, with a looser criteria than that stated above, which 
are constrained to come from a single vertex. We also require that the invariant mass be 
within 5 times the width of the A mass peak, which has an r.m.s. width of 1.4 MeV. 

A. Efficiency Determinations 

To derive event selection efficiencies we simulated hadronic events using the JETSET 7.3 
qq event generator 0, then followed through the full GEANT 3.21-based Q CLEO-III 
detector simulation. For five-flavor hadronic and A°A° events in the A& scan region, we 
generated Monte Carlo samples using the same generator with the properties described in 
section |H] The efficiencies obtained from these simulations are presented in Table |Hj where 
we list the both the hadronic event selection efficiency and the efficiency for detecting a 
hadronic event with an anti-proton. These efficiencies include the branching ratios for the 
various processes into anti-protons in the second column. We take £>(A£ — > pX) = 0.50. The 
row for bb includes only B meson production with additional pions allowed. As one would 
expect, the efficiencies for bb and A°A° are very similar. The slightly lower efficiency for 
A°A° arises from higher average jettiness for A°A° events. 

TABLE II: Selection efficiencies for hadronic events and those with anti-protons. 
Data samples Selection efficiency for Selection efficiency for 



hadronic events (%) hadronic events with an p (%) 




89.9 ± 1.2 ± 2.7 
28.1 ± 2.5 ± 0.8 
0.024 ± 0.005 ± 0.001 



25.5 ± 0.2 ± 0.8 2.1 ± 0.1 ± 0.1 

85.5 ± 0.9 ± 2.6 26.8 ± 0.1 ± 5.4 

21.9 ± 0.4 ± 0.7 1.8 ± 0.2 ± 0.1 



4.0 ± 0.2 ± 0.3 
2.0 ± 0.3 ± 0.2 

< 10~ 5 



TT 
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The errors listed in Table ITU are statistical and systematic, respectively. The systematic 
error for the hadronic event selection requirement is estimated from the variation in the 
number of hadronic events (corrected by efficiency and background) when changing selection 
requirements. The systematic error for the proton identification has been evaluated from 
proton efficiency measurements using reconstructed A events from data and then comparing 
with the equivalent MC estimation. 

Our simulations also give us the selection efficiency for detecting an event containing 
either a A or an A from A^A^ decay of 16.6±0.llJ;o%, including the B(A — ► pn~). Note that 
the PDG world average for B(A C — > p anything) is (50 ±16)%. Similarly B(A C — > A anything) 
is (35 ± 11)% The errors on these rates will be included separately as systematic effects. 

B. Systematic Errors 

The systematic errors in determining A°A° production are given in Table ITTT1 The largest 
error is due to the unknown branching fraction of B(A C — > pX) to which we assign a 32% 
error. We also include errors on the hadron selection efficiency and the background in 
the hadronic event sample, evaluated by varying our selection criteria as well as taking 
into account the variation with Ecm, the anti-proton identification efficiency evaluated by 
examining a larger sample of A — > pn~ data, and the luminosity measurement uncertainty 
estimated as 1% 

The total systematic error found by adding these elements in quadrature is 2.7%, 32% and 
31% on the determination of A\A\ production using bb, anti-protons and A's, respectively. 

TABLE III: List of systematic errors in determining AjjAjJ production 



Source Error 

Hadron efficiency ±3 

A£ — ► A+A branching ratio ±4 

Proton identification efficiency ±4 

A^~ — > pX branching fraction ±32 

A+ — > A A branching fraction ±31 
Total background of hadronic events ±2 

Luminosity ±1 



IV. THE ESTIMATED bb CROSS SECTION 

The hadronic cross section is generally expressed in terms of its ratio R to the point cross 
section e + e~ —>■ p^pT. To search for resonant or non- resonant production of A\A\ in e + e~ 
collisions we measure the bb cross section over the energy range of the scan. Theoretically, 
R b ^ can be expressed as follows: 



R bb ~ R bE 



l + a s /n + C 2 (a s /7i) 2 + C 3 (a s /7cf , (2) 
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where = N c qb. N c is the number of quark colors, q b is the b quark charge and a s is the 
strong coupling constant. The constants are C2 = 1.409 and C3 = —11.767 [14] . In our 
energy regime, we expect a value for R bb of 0.35. 

To find the bb cross section we subtract the R2 four-flavor continuum data distribution 
from the higher energy data, correct for the efficiency of the R2 cut and the hadronic selection 
criteria and divide by the relevant luminosity. We use a value of the cross-section for e + e~ — > 
equal to 86.8 nb/s, where s is the square of the center-of-mass energy in units of GeV. 
However, we do not make a precise measurement of bb cross section due to uncertainties 
in the correct scaling factors of two-photon events and initial state radiation contributions 
in different energy regions. Here we wish to measure any possible enhancement above the 
A°A° threshold. Our results are presented in Fig. 0(a). 



0.36 



0.32 



EC 0.28 



0.24 



0.20 



a) 



V 



v v 



_L 



_L 



0.36 



0.32 



cc 0.28 



0.24 



0.20 




11.15 11.20 11.25 11.30 11.35 

s 1/2 (GeV) 



11.15 11.20 11.25 11.30 11.35 

s 1/2 (GeV) 



FIG. 5: The estimated bb cross section in units of R. The error bars on the data points represent 
both the statistical and the systematic errors summed in quadrature, (a) The solid line shows a fit 
to a horizontal line, (b) The solid line shows a fit to Equation |3J The fits are described in the text. 



V. UPPER LIMITS ON A b PRODUCTION 

In this energy regime we expect that the R value will be constant in the absence of any 
resonant or threshold increase due to A°A° production. There are no statistically significant 
excesses above a constant value of R, suggesting no resonant production of bb types of events. 
There is an important caveat concerning the limit using the bb cross-section. It may very 
well be that opening up the A°A° channel comes at the expense of a lower in rate of other 
channels so that the total bb rate remains constant. Should this occur our limit, in this (bb ) 
case, would be meaningless. In fact, a fit to flat line for bb yields a \ 2 °f 14.2 for 29 degrees 
of freedom. This fit is shown on Fig. |S] (a). 

We can look for an increase in A°A° production that mimics the threshold turn on as 
a function of center-of-mass energy of e + e~ — > t + t~. The line in Fig. |5] (b) represents a 
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two- component fit. The first component is a straight line without any slope allowed up to a 
Eqm of 11.24 GeV, twice the A& mass. The second component uses a shape similar to one 
proposed by the BES collaboration [l5|, but simplified by explicitly calculating the Coulomb 
interaction and final state radiation; the final form of this function is: 

a{s) = Ax6{^-2m{Al)){^-2m{A Q b )) om + R , (3) 

where A is a fit parameter, 9{y) is step function, for y<0 and 1 for y> 0, m(A b ) is the mass 
and Rq is the observed cross section below threshold. (We are assuming this form applies 
only near threshold.) 



1631004-101 
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FIG. 6: The cross section for events with at least one anti-proton normalized by <r(e + e~ — > 

(The data have not been corrected for hadronic event efficiencies.) The solid lines show fits to 

Equation |3J The errors are statistical only. 

The cross sections for events with anti-protons are shown in Fig. |U1 The data have been 
corrected for the momentum dependent efficiency of identifying anti-protons, but not for 
hadronic event selection. The data are fit to the BES function given in Eq. ©• The fitted 
parameters used to set upper limits are listed in Table II VI 

The cross sections for events with A's are shown in Fig. [7| The data have been corrected 
for both A reconstruction efficiency and the branching ratio for A b A b into A plus A. The 
fit to data uses the BES function given in Eq. (jSJ). The fitted parameters used to set upper 
limits are listed in Table HVl 

There is no significant resonance peak in the scan range, nor any evidence for a growth 
above threshold. Using these fits we calculate 95% confidence level upper limits for 
A b A b production above threshold, as shown in Fig. |HJ Here we take the upper limit as 

a(s)r er = (A + 1.64 x SAi) * (^i - 2m(A°)) 062 /e, , (4) 

where Aj is the fit value from Table HV| 5Ai is its error, is the relative A^ efficiency for each 
of the three different methods of 0.95, 0.29, and 0.86, for bb, p and A searches, respectively. 
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0.1 OU i i i i I i i i i I i i i i I i i i i I i i l 

11.15 11.20 11.25 11.30 11.35 

s 1/2 (GeV) 

FIG. 7: The cross section for events with at least one A normalized by a{e + e~ — ► (The 
data have not been corrected for hadronic event efficiencies.) The solid lines show fits to Equation^ 
The errors are statistical only. 

The 0.95 results from the relative efficiency of continuum bb production to A&A;,, the 0.27 
is the product of the A^Afe decay rate into anti-protons and the efficiency of the hadronic 
event selection, and the 0.86 is hadronic event selection for A b A&. The systematic errors 
are included only in the limits using bb production. In the other two cases the systematic 
errors on the inclusive p and A branching ratios worsen the upper limits by 32% and 31%, 
respectively. 

We determine upper limits for production of a resonance that would decay into A°A° , 
similar in spirit to T(4S) — > BB. Here we take two possible intervals for either a narrow 6 
MeV wide resonance or a wider, arbitrarily chosen, 18 MeV resonance. For the first case we 
fit a horizontal line to our data up to the A°A° threshold of 11.24 GeV and then estimate 
the upper limit for a cross-section excess in each 6 MeV interval of center of mass energy. 
These 95% confidence level upper limits are shown in Fig. EJa). 

For the second case, we fit all our data to a horizontal straight line while excluding an 
18 MeV wide interval of center-of-mass energy. We then calculate the 95% confidence level 
upper limit by calculating the difference of the data relative to the fit line. These limits are 



TABLE IV: Numerical values of parameters found by fitting Eq. El to our data. 



Selection criteria 


At 




Rot 


bb 


(0.21 ±3.82) 


x 10" 2 


0.322 ± 0.007 


Anti-proton 


(0.84 ± 1.20) 


x 10~ 2 


0.333 ± 0.002 


A 


(0.15 ±5.49) 


x 10" 2 


0.201 ± 0.010 



Twice the A& mass is fixed to 11.24 GeV. 
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shown in Fig. E^b). 




FIG. 8: The fractional upper limits at 95% c. 1. for A^A^ production obtained using A (solid line), 
anti-proton (dashed line) and the bb (dotted line) yields set by using the BES function. For the 
bb case only, systematic errors have been included. 

No resonant enhancement reminiscent of the T(4S) resonance is observed. Using the 
threshold function we can set an upper limit at our highest energy point of 11.383 GeV 
on the ratio of A° to bb production. These limits are given in Table For bb production 
we use two values - the first is as defined in Eq. [21 the second is determined by fitting 
R bb values assuming no enhancement along scan range. These values are = 1/3 and 
R bl = 0.322 ± 0.004. 

The limits based on this function become lower toward lower energy as we approach the 
production threshold. The anti-proton and A samples are somewhat correlated in that anti- 
protons from A decay are often included in both samples, so we choose not to combine these 
limits. 

TABLE V: Upper limits at 95% c.l. on the ratio of A£A£ to bb production at 11.383 GeV. 



method 


95% c.l. (statistical only) 95% c.l. (statistical & systematic) 




R A°I° b / R bb R A°I° b /Rbl R A°X /R bb 


R A°A°J Rb ~ b 


bb 


6.0% 


6.2% 


Anti-proton 


9.2% 9.5% 12.2% 


12.5% 


A 


9.9% 10.2% 12.9% 


13.3% 



VI. CONCLUSIONS 

We do not observe any resonant or threshold enhancement of A b A b production in the 
center of mass energy region just above threshold, resulting in 95% confidence level upper 
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11.15 11.20 11.25 11.30 11.35 11.21 11.25 11.29 11.33 11.37 

s 1/2 (GeV) S 1/2 (GeV) 

FIG. 9: Upper limits at 95% c. 1. for A°Ajj production obtained using A (solid line), anti-proton 
(dashed line) and the bb (dotted line) yields, (a) The upper limits have been set in 6 MeV center 
of mass energy intervals in the scan region, (b) Upper limits in 18 MeV wide intervals. For the 
bb case only, systematic errors have been included. 



limits on the order of 0.05-0.10 units of R. The 95% confidence level upper limits from anti- 
proton and A production are 12.8% and 12.9% of R®r, respectively, at our highest energy 
point if they are modelled as a growth above threshold. In order to effectively study A5 
decays at e + e~ machines, it may be necessary to go to higher center of mass energies. 
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